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Structural Transition in Layered Asi_ x P x Compounds: A Computational Study 
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As a way to further improve the electronic properties of group V layered semiconductors, we 
propose to form in-layer 2D heterostructures of black phosphorus and grey arsenic. We use ab initio 
density functional theory to optimize the geometry, determine the electronic structure, and identify 
the most stable allotropes as a function of composition. Since pure black phosphorus and pure grey 
arsenic monolayers differ in their equilibrium structure, we predict a structural transition and a 
change in frontier states, including a change from a direct-gap to an indirect-gap semiconductor, 
with changing composition. 

PACS numbers: 73.20.At, 73.61.Cw, 61.46.-w, 73.22.-f 


INTRODUCTION 

Few-layer structures of group V elements, includ¬ 
ing phosphorene and arsenene, are emerging as promis¬ 
ing candidates for two-dimensional (2D) electronic ma¬ 
terials application [JH5]. Different from semi-metallic 
graphene [7], these systems display a large band gap 
while still maintaining a high carrier mobility [3j I8FHT] . 
Even though phosphorus and arsenic are both group V 
elements, they crystallize in different structures. Most 
stable and thus more abundant are the layered allotropes 
such as black phosphorus [12], with the designation A17 
or <a-P, and grey arsenic [4], with the designation A7 or 
/3-As. [13] Since a conversion of the a to the /3 phase is 
possible [HI [15], combining both elements in the same 
layer and changing the composition is bound to cause 
a structural transition[16, Uj. Since both structures 
are almost equally stable, we may expect phase coex¬ 
istence that should bring an unexpected richness in both 
structural and electronic properties [5, 16, [18] 191. This 
way of isoelectronic doping could turn into an effective 
way to fine tune the electronic properties, improve the 
carrier mobility [5] and reduce the chemical reactivity of 
the compound from those of pristine of phosphorene and 
arsenene ]20]. 

In this study we report ab initio density functional the¬ 
ory (DFT) calculations of Asi-^P^ monolayers. We de¬ 
termine the optimum geometry, the electronic structure, 
and identify the most stable allotropes as a function of 
composition. We predict a structural transition from the 
a to the /3 phase to occur near x=0.93. This structural 
transition is accompanied by a change in frontier states 
from lone pair electron states in a-P to cr-bond states in 
yd- As and from a direct to an indirect fundamental gap. 

RESULTS AND DISCUSSION 

We present in the following computational results for 
the equilibrium geometry and electronic structure of 
Asi_*P x compounds as a function of composition. The 
monolayer structures have been optimized using DFT 


with the Perdew-Burke-Ernzerhof (PBE)[2T] exchange- 
correlation functional, as discussed in the Methods Sec¬ 
tion. Since group V elements phosphorus and arsenic 
are both threefold coordinated and display a tetrahedral 
bonding character, they possess the freedom to arrange 
atoms in a layered structure that is not flat, which leads 
to various buckled allotropes that are topologically re¬ 
lated to a honeycomb lattice[22]. Among these, the a 
and yd phases are most stable [14]. Under ambient pres¬ 
sure, phosphorus favors the a and arsenic the (3 struc¬ 
ture, but an <a—>yd transition has been reported under 
high pressure[fl5l [23] . As seen in Fig. [lja), a monolayer 
of a-P (black phosphorus) can be viewed as a distorted 
honeycomb lattice with a rectangular primitive unit cell. 
The structure of a yd- As (grey arsenic) monolayer, pre¬ 
sented in Fig. 0b), resembles more closely the honey¬ 
comb lattice of graphene with a hexagonal unit cell. 

We explored the structures of Asi-^P^ compounds 
with 0 < x < 1 by considering different arrangements of 
As and P atoms in systems with x = 0.25, 0.5 and 0.75. 
We optimized the lattice structure for each system and 
found that the lattice constants depend primarily on x 
and increase with increasing As concentration due to the 
larger atomic radius of arsenic in comparison to phos¬ 
phorus. In the strongly anisotropic a phase, we found 
an increase of 4% along the softer ai direction and a 
10% increase along the harder <22 direction when moving 
from pristine phosphorus to arsenic. The corresponding 
change in both a\ and a 2 was an increase by 9% in the 
isotropic /3 phase. The data for the dependence of the 
lattice parameters on composition are presented in the 
Supporting Information. 

Since the cohesive energy of phosphorus is larger than 
that of arsenic, we expect a larger cohesive energy in 
phosphorus-rich structures. In the first approximation, 
ignoring structural and short-range order differences, we 
expect the cohesive energy of the Asi-^P^ compound 
to be a linear combination of the cohesive energies of 
the pristine components in their respective structures, 
^ P3 C oh — (1 •E)P3 C oh ( As)-\-xE co ji (P)• In n com¬ 
pound with a particular structure, we define the co¬ 
hesive energy per “average atom” by E co h(Asi- x P x )= 
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FIG. 1. (Color online) Equilibrium structure of a monolayer of (a) black phosphorus (a-P), (b) grey arsenic ((3- As), and (c) the 
/3 allotrope of the Aso.25Po.75 compound. Ball-and-stick models in top and side view, with the primitive unit cell highlighted 
by shading, are shown in the upper panels and the Brillouin zones in the lower panels. 


-E tot (As 1 - x P x )/N + (1 - x)E tot ( As akom)-\-xE tot (P 
atom), where E tot (Asi- x P x ) is the total energy of the 
AT-atom unit cell, E tot ( As atom) is the total energy of 
an isolated As atom and E to t( P atom) the total energy 
of a P atom. As seen in Fig. [2ja) , the cohesive energy 
E co h{ Asi-^P^) will generally deviate from the expecta¬ 
tion value < E co h > for specific geometries and atomic 
arrangements, but the deviations are rather small. 

To better investigate these deviations, we de¬ 
fine the relative stability of a given structure by 
E C oh (Asi—^P x) — < E co h > and show the results in 
Fig. [2jb) . In general, we find that combining both el¬ 
ements in the same layer is always associated with an en¬ 
ergy penalty with respect to the < E co h > value, caused 
by the size mismatch of As and P atoms. For a given 
composition, we found cohesive energy differences be¬ 
tween particular arrangements of As and P atoms to be 
as large as 20 meV/atom at x = 0.5. Since the cohe¬ 
sive energy ranges of the different atomic arrangements 
within the phases overlap, we expect a coexistence of the 
two phases in realistic samples. 

We found the monolayer of /3-As to be more stable by 
~26 meV/atom than that of a-As, whereas the mono- 
layer of a-P is favored by ^3 meV/atom over (3- P. As 
seen in Fig. [2jb) , the energy penalty in the /? phase is 
smaller than that in the a phase for lower x values. Only 
in very phosphorus-rich compounds we find the a phase 
to be more stable. Consequently, we expect a structural 
transition between the a and the (3 phase to occur as a 


function of composition. Our results indicate that this 
structural transition should occur near x~0.93. 

To illustrate how the relative position of arsenic and 
phosphorus atoms affects the stability of the compounds, 
we present three different structures of AS 0 . 5 P 0.5 in the a 
phase in Fig. [2j(c) to (e)] and in the (3 phase in Fig. [2j(f ) 
to (h)]. As the cohesive energy of phosphorus is larger 
than that of arsenic, P-P bonds are stronger than As- 
As bonds. Consequently, the structure with the largest 
number of P-P bonds should be most stable. On the 
other hand, maximizing the P-P interaction would im¬ 
ply segregating P from As, which will strain the struc¬ 
ture. There is a trade-off between the two trends. For 
the sake of convenience, we called the least stable struc¬ 
tural arrangement “structure I” and the most stable ar¬ 
rangement “structure III” in Fig. |2j As seen in Fig. [2je) 
and[2jh), the two most stable structures in either phase 
contain an alternating arrangement of isolated P-P and 
As-As dimers. In this arrangement, the strain caused by 
different bond lengths can be minimized, while still keep¬ 
ing at least some P-P bonds. Strain is also low in the a-II 
and (3-1 structures, but absence of P-P bonds makes them 
less stable. The largest number of P-P bonds is realized 
in the a-1 and (3-11 structures, but the large strain en¬ 
ergy caused by the coexistence of phosphorus and arsenic 
chain arrangements makes these structures energetically 
unfavorable. 

As mentioned above, the relatively small stability dif¬ 
ference between the a and (3 phase of Asi-^Pa, com- 
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FIG. 2. (Color online) (a) Cohesive energy E co h of an Asi-ccPcc monolayer as a function of the phosphorus concentration x. 
The most stable allotropes are shown by larger symbols. < E co h >, shown by the dash-dotted line, is a linear interpolation 
between the most stable allotropes of As and P. (b) Relative stability E co h~ < E co h > of Asi-ccPcc as a function of x. The 
vertical dashed line indicates the expected composition for an a^/3 structural transition. Ball-and-stick models of selected 
stable structures of AS0.5P0.5 in the [(c) to (e)] a phase and [(f) to (h)] (3 phase. For each phase, the structural indices I, II 
and III are arranged in the order of increasing stability. 
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FIG. 3. (Color online) Fundamental electronic band gap E g in Asi- x Pcc monolayers. DFT-PBE values of E g are presented as 
a function of composition for compounds in the (a) a- and (b) /3-phase. The dashed lines are guides to the eye. 


pounds suggests a likely coexistence of different al¬ 
lotropes. Understanding the electronic properties of such 
a complex system requires obtaining information about 
every such structural arrangement. We have performed 
the corresponding calculations and show the fundamental 
electronic band gap of Asi-^P^ compounds as a function 


of composition in Fig. [3j As indicated in Fig. [3][a), the 
electronic band gap value of the a phase does not de¬ 
pend sensitively on the phosphorus concentration x and 
lies in the range between 0.9 — 1.0 eV. For the /3 phase of 
Asi-^P^ compounds, the band gap is much larger and 
it’s value increases with increasing phosphorus concentra- 
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FIG. 4. (Color online) Electronic band structure of Asi-^Pcc monolayers in the at-phase [(a) to (e)] and the /3-phase [(f) to (j)]. 
The structural arrangements for each composition are shown in the insets. The changing role of the valence band maxima, 
shown by the dashed circles labeled “A” at T and “B” off-T, is discussed in the text. The position of the conduction band 
minima is indicated by the dashed circles labeled “C”. 


tion, as seen in Fig. ^b). The band gap values not only 
display a larger value range from 1.6 — 2.0 eV, but also 
depend on the relative arrangement of P and As atoms 
within the unit cell at a given composition. 

More interesting than the absolute value of the band 
gap are changes in the band structure of Asi-^P^ com¬ 
pounds caused by changing composition, which should 
be correctly captured by DFT-PBE. Especially impor¬ 
tant appears to be the position of the valence band max¬ 
imum, since these compounds are expected to behave as 
p-type semiconductors. The band structure of ot-phase 
compounds is shown in Fig.[4][a) for pristine phosphorus, 
Fig. Be) for pristine arsenic, and in Figs, ^b) through 
^d) for intermediate compositions. Careful comparison 
of these band structure results indicates a transition from 
a direct gap in a-P to an indirect gap in a- As. For the 
sake of discussion, we found it useful to identify the va¬ 
lence band maximum at T as point “A”, another local va¬ 
lence band maximum along the T — X line as point “B”, 
and the conduction band minimum at T as point “C”. 
With increasing concentration of As, the valence band 
maximum switches from “A” to “B”, with the transition 
occurring near Aso. 75 Po .25 as seen in Fig. ^d). Con¬ 
sequently, the character of the fundamental band gap in 
Asi-^P^ compounds in the a phase is expected to change 


from direct for x>0.75 to indirect for x<0.75. This be¬ 
havior is reminiscent of the direct-to-indirect gap transi¬ 
tion in a pure <a-phosphorene monolayer that is induced 
by tensile in-layer strain and may be rationalized by the 
fact that the larger atomic radius of As atoms causes such 
strain. In reality, we find that the relative positions of the 
“A” and “B” peaks depend not only on the composition, 
but also the relative arrangement of phosphorus and ar¬ 
senic atoms. We compare the electronic band structure 
of different structural arrangements with x=0.5 in the 
Supporting Information. 

The band structure changes in /3-Asi- x P x compounds 
with changing composition are presented in Fig. |4][f) to 
Bit Similar to the a-phase structures discussed above, 
we label the valence band maximum at T as point “A”, 
another local valence band maximum away from T as 
point “B”, and the bottom of the conduction band as 
point “C”. Even though - unlike the <a-phase compounds 
- all /3-phase structures are indirect-gap semiconductors, 
the position of the valence band maxima and conduction 
band minima still changes with composition. There is a 
change in /3-Asi- x P x from the top of the valence band 
being at “B” in P-rich compounds with £>0.5 to “A” at T 
in As-rich compounds with £<0.5. Similar changes in the 
valence frontier states have been observed in a pure As 
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monolayer under strain^]. Also here, we find that point 
“A” corresponds to in-plane cr-bond states, whereas point 
“B” corresponds to lone pair electron states, [4] same as 
in a /3-P monolayer. In other words, changing the com¬ 
position allows to change the character of the frontier 
states in the valence band region and thus to effectively 
tune the electronic properties of the system. 

In conclusion, we have performed ab initio density 
functional theory (DFT) calculations of Asi_ x P x mono- 
layers in search of ways to further improve the electronic 
properties of group V layered semiconductors and reduce 
their chemical reactivity. We have determined the opti¬ 
mum geometry, the electronic structure, and have iden¬ 
tified the most stable allotropes as a function of compo¬ 
sition. We have found the most stable allotropes to be 
based on the a (A17) structure of black phosphorus and 
the (3 (A7) structure of grey arsenic. Since the stabil¬ 
ity difference between these two phases is very small, we 
expect coexistence of a- and /3-type structures within a 
large composition range. Based on our results, we ex¬ 
pect a structural transition from the a to the (3 phase to 
occur near x=0.93. This structural transition should be 
accompanied by a change in frontier states from lone-pair 
electron states in a-P to cr-bond states in (3 -As and from 
a direct to an indirect fundamental gap. 


METHODS 

Our computational approach to gain insight into the 
equilibrium structure, stability and electronic proper¬ 
ties of various phosphorene structures is based on ab 
initio density functional theory as implemented in the 
SIESTA code.|24 : We used periodic boundary condi¬ 
tions throughout the study. We used the Perdew-Burke- 
Ernzerhof|2T] exchange-correlation functional, norm- 
conserving Troullier-Martins pseudopotentials [25] . and a 
double-^ basis including polarization orbitals. The recip¬ 
rocal space was sampled by a fine grid[26 of 8x8x1 k- 
points in the Brillouin zone of the primitive unit cell. We 
used a mesh cutoff energy of 180 Ry to determine the 
self-consistent charge density, which provided us with a 
precision in total energy of < 2 meV/atom. All geome¬ 
tries have been optimized by SIESTA using the con¬ 
jugate gradient method[27], until none of the residual 
Hellmann-Feynman forces exceeded 10 -2 eV/A. 

We thank Bilu Liu for useful discussions. This study 
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thors declare no competing financial interest. 
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